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ABSTRACT 

We analyze published reverberation mapping data for three Seyfert galaxies (NGC 3227, 
NGC 3516, and NGC 4593) to refine the mass estimate for the supermassive black hole in the 
center of each object. Treatment of the data in a manner more consistent with other large 
compilations of such masses allows us to more securely compare our results to wider samples of 
data, e.g., in the investigation of the M^-cr* relationship for active and quiescent galaxies. 



Subject headings: galaxies: active - 
nuclei — galaxies: Seyfert 

INTRODUCTION 



galaxies: individual (NGC 3227, NGC 3516, NGC 4593) — galaxies: 



Reverberation mapping (Blandford & McKee 
1982; Peterson 1993) is an important technique for 
probing the supermassive black holes (SMBHs) at 
the centers of active galactic nuclei (AGNs) . Mea- 
suring the time lags, t, between fluctuations in 
the highly variable AGN continuum and the vari- 
ations of emission lines arising from the broad- line 
region (BLR) provides a responsivity-weighted ra- 
dius for the BLR gas distribution around the 
SMBH. When emission-line time lags are com- 
bined with calculations of the gas velocity width, 
cr, under the assumption of virialized gas motions, 
it becomes possible to estimate the SMBH mass. 

Multiple lines of evidence suggest that to as- 
sume a virialized BLR is a reasonable approx- 
imation. In AGNs for which several emission 



1 Department of Astronomy, The Ohio State 
University, 140 West 18th Avenue, Colum- 
bus, OH 43210; onken@astronomy.ohio-state.edu, 
peterson@astronomy.ohio-state.edu 

2 Department of Astronomy, University of Florida, 211 
Bryant Space Science Center, Gainesville, FL, 32611-2055. 
Current address: Department of Physics and Astronomy, 
Georgia State University, One Park Place South SE, At- 
lanta, GA, 30303; dietrich@chara.gsu.edu 

3 Division of Physics and Astronomy, Department of 
Physical Sciences, University of Hertfordshire, College 
Lane, Hatfield, Herts ALIO 9AB, United Kingdom; 
ar@star.herts.ac.uk 

4 Anton Pannckock Institute, University of Amsterdam, 
Kruislaan 403, Amsterdam 1098 SJ, The Netherlands; is- 
abel@science.uva. nl 



lines have been reverberation-mapped (NGC 5548, 
3C 390.3, NGC 7469, NGC 3783), an inverse re- 
lationship is found between r and a that is con- 
sistent with expectations from virialization (Pe- 
terson & Wandel 1999, 2000; Onken & Peterson 
2002). In addition, the close concordance between 
AGNs and quiescent galaxies in the M^-a* plane 
under the virial assumption (relating black hole 
mass to the stellar velocity dispersion of the galac- 
tic spheroid; Ferrarese et al. 2001) would seem to 
imply that any systematic error in reverberation 
masses is small (a factor of ^3 or less). Krolik 
(2001) estimated a similar level of expected sys- 
tematic scatter assuming gravitational domination 
of the BLR. 

Most of the statistical studies of reverberation- 
based black-hole masses use the homogeneous 
compilations of Wandel, Peterson, & Malkan 
(1999, hereafter WPM) and Kaspi et al. (2000). 
Missing from these compilations are AGNs ob- 
served by the "Lovers of Active Galaxies (LAG)" 
collaboration that monitored several AGNs in 
early 1990 (Robinson 1994). Three of the LAG 
sources, NGC 3227 (Salamanca et al. 1994), 
NGC 3516 (Wanders et al. 1993), and NGC 4593 
(Dietrich et al. 1994), have well-determined 
emission-line lags and certainly meet the quality 
criteria for inclusion in these studies, and indeed 
were included in the compilation of Ho (1999; 
based on H/3 only). These galaxies were omit- 
ted from the WPM compilation only because the 
data had not been analyzed in the same fash- 
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ion as the other sources; specifically, (a) WPM 
use the model-independent Monte Carlo method 
of Peterson et al. (1998) to assess uncertain- 
ties in emission-line lags, and (b) WPM use the 
FWHM of the emission-line in the root-mean- 
square (rms) spectrum, rather than the mean 
spectrum, to characterize the BLR line-of-sight 
velocity width, AV, that is used to form the virial 
product (AV) 2 t. Because we have converted both 
our time lags and velocity widths to the AGN rest 
frame, our results supplement the compilation of 
Kaspi et al. (2000). (WPM did not correct for the 
time dilation of their r values.) 

In this contribution, we reanalyze the LAG 
data on these three sources in the same fashion 
as Kaspi et al. (2000), with the intent of enlarging 
that homogeneous database. The total number of 
AGNs which have been reverberation-mapped (a 
few dozen) is small enough that additional objects 
will assist in on-going statistical investigations of 
(what are hoped to be) fundamental physical re- 
lationships. The main aspect on which we have 
worked is the careful construction the rms spectra 
over the spans of the respective observing cam- 
paigns. With the rms spectra for these AGNs and 
updated time lag values and uncertainties, we are 
able to derive masses for the SMBHs in these three 
galaxies. 

In § 2, we briefly describe the original observing 
campaigns and our data reduction. We explain the 
details of our time series and velocity width anal- 
yses in § 3. Our estimates of the SMBH mass in 
each galaxy are presented in § 4. In § 5, we address 
the Mbh-cr* relationship, and we then summarize 
our conclusions (§6). 

2. OBSERVATIONS AND DATA RE- 
DUCTION 

The LAG consortium examined the Ha and H/3 
wavelength regimes with long slit spectroscopy of 
moderate resolution and with the temporal sam- 
pling rates listed in Table 1. The line marked 
"cont." in Table 1 for NGC 3516 indicates the 
sampling statistics derived for the combination of 
imaging and spectroscopic measurements of the 
~4000-5000A continuum. The continuum statis- 
tics for NGC 3227 and NGC 4593 arc identical to 
the H/3 and Ha entries, respectively. F var indi- 
cates the "excess variance," the mean fractional 



variation of each dataset (see Rodrfguez-Pascual 
et al. 1997). The ratio of the maximum to mini- 
mum flux value is given in Table 1 as R max - 

In general, no alterations were applied to the 
published light curve data. However, in order to 
investigate the velocity width characteristics of the 
broad emission lines, we used the wavelength re- 
gion around the narrow lines ([O III] and [S II], 
which should only vary on timescales much longer 
than the campaigns) to scale the spectra (van 
Groningen & Wanders 1992). The absolute scal- 
ing was done so as to yield mean narrow-line fluxes 
equal to the values derived in the original papers. 
The slit width selected was an attempt to balance 
spectral resolution and total detected flux. How- 
ever, the small slit width introduced a non-trivial 
degree of difficulty in the removal of seeing effects, 
as noted below. 

2.1. NGC 3227 

The LAG campaign on NGC 3227 was de- 
scribed by Salamanca et al. (1994). The analysis 
was complicated by the fact that the two meth- 
ods they employed for removing the contribution 
of the stellar host gave significantly different an- 
swers: they derived a 24% contribution to the 
continuum by examining the magnesium triplet of 
a template galaxy of the same Hubble type, and 
a 40% continuum contribution from comparisons 
with the magnesium triplet absorption of the bulge 
of NGC 3227. We analyzed the data with each of 
the potential stellar levels removed and found the 
differences to be insignificant. 

A detailed investigation of various origins of un- 
certainty in their flux calibration was conducted 
by a portion of the LAG group (Baribaud et al. 
1994). In addition to the effects of slit miscenter- 
ing and seeing, they estimate the uncertainties in- 
curred by their scaling of the spectra by the narrow 
lines, and conclude an overall flux error of 5-8%. 
We have assumed the more conservative 8% errors 
for all of the flux data. Our analysis made use of 
the full LAG dataset, consisting of 23 Ha obser- 
vations and 14 H/3 observations. The light curve 
cross-correlations were made with respect to the 
continuum at 5020A. 

The systemic redshift of NGC 3227 was taken 
to be z = 0.00380 ± 0.00002 (Keel 1996). 
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2.2. NGC 3516 

Spectra of NGC 3516 were taken at 22 separate 
epochs for Ha and 21 epochs for H/3. The con- 
tinuum light curve for this object was constructed 
from a combination of B-band imaging and spec- 
tral fits to the power-law continuum at 4905A. 
The results of the LAG study have been published 
by Wanders et al. (1993) and Wanders & Home 
(1994). The spectra we obtained for this object 
had already been corrected for seeing, the stellar 
flux from the host galaxy, and contributions from 
narrow lines (see Wanders et al. 1992). In addi- 
tion, the [S II] and [O III] narrow lines had been 
subtracted from the Ha and H/3 spectra, respec- 
tively, preventing any rescaling of the flux levels. 
Several observations were removed from the calcu- 
lation of the rms spectrum for a variety of reasons 
(extremely poor seeing, object miscentering, and 
other anomalies). Among the Ha dataset, we have 
excised observations from the following dates: 31 
January, 16 February, 15 April, and 13 May. From 
the H/3 analysis we have removed 24 March, 15 
April, and 13 May. 

Conversion of time lags and line-of-sight veloc- 
ity widths to the rest frame was made with z = 
0.00884 ± 0.00002 (Keel 1996). 

2.3. NGC 4593 

Spectra of NGC 4593 were obtained at 22 
epochs around Ha and at 11 epochs around H/3, 
and the continuum used was at 6310A. The Ha 
line was contaminated by strong [N II] AA6548, 
6584 emission, as well as a narrow Ha contribu- 
tion. Using the shape of the narrow [O III] A5007 
line as a template for all of the narrow lines, and 
the flux ratios between the [N II] lines and narrow 
Ha that Dietrich et al. (1994) calculated, we have 
bracketed the amount of narrow line removal by 
limiting cases in which too much or too little nar- 
row line flux was subtracted. Conducting this for 
each Ha spectrum, we then separately analyzed 
our oversubtracted and undersubtracted datasets. 
The AU difference between the oversubtraction 
and undersubtraction was on the order of 1%, 
much smaller than the AV uncertainties as well 
as the range of values derived by doing no rescal- 
ing of the spectra or no removal of the narrow-line 
contribution. 

We converted our results to the galaxy's rest 



frame using z = 0.00900 ± 0.00013 (Strauss et al. 
1992). 

3. DATA ANALYSIS 

Using the published light curves, we generated 
the cross-correlation functions (CCFs) between 
the continua and the emission lines with the inter- 
polated CCF technique as implemented by White 
& Peterson (1994). We calculated both the peak 
of each CCF (T pea k) and its centroid (r cent ; above 
a threshold of 80% of the peak value correlation 
coefficient) with an interpolation unit of 0.1 days. 

Uncertainties in the time lag values were esti- 
mated with the method described by Peterson et 
al. (1998). This model-independent technique ad- 
dresses two major sources of uncertainty in the 
CCFs. First, the issue of flux errors is handled by 
adding to each flux value a randomly generated 
deviation. This deviation is computed such that a 
large number of realizations will produce a Gaus- 
sian scatter with a standard deviation equal to 
the quoted uncertainty (this procedure is dubbed 
"flux randomization (FR)"). Second, the poten- 
tial errors arising from the limited temporal sam- 
pling of the observational campaign are accounted 
for by "random subset selection (RSS)." The RSS 
method consists of randomly selecting N data 
points from a sample of size N, where each se- 
lection draws from the entire sample and repeated 
selections are ignored. This allows us to "boot- 
strap" our way to reasonable error estimates, while 
preserving the underlying form of the time series. 
The FR/RSS method was implemented 1000 times 
for each continuum-line combination to build up a 
distribution of the peak and centroid time lags, de- 
noted generally as a cross-correlation peak distri- 
bution (CCPD; Maoz & Netzer 1989). We use the 
CCPD to define the time lag values within which 
68% of the Monte Carlo results lie, and these val- 
ues set our time lag uncertainties. Table 2 and 
Figure 3 show our cross-correlation results. 

The procedure we applied to the mean and rms 
spectra of each AGN was to make two estimates 
of the continuum, at the extrema of where the 
level could potentially be set. The two resulting 
measurements of the velocity full-width at half- 
maximum, AV, were then averaged and the dif- 
ference taken to be twice the one-sigma errors. 
We give the results of our AV measurements from 
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both the mean and rms spectra in Table 3. 

4. MASS DERIVATIONS 

For consistency with previously published com- 
pilations, we calculate the reverberation mass as 



M r 



3 c t (AV) 2 
4 G 



(1) 



where AV is taken from the rms spectrum, c is 
the speed of light, and G is the gravitational con- 
stant. The asymmetric errors in the time lags are 
propagated through to the masses derived from 
the individual lines. Then, the masses for each 
line are combined in a weighted average. We have 
measured the reverberation masses for our three 
targets and present the results in Table 4. Object- 
specific comments and comparisons with previous 
results follow. 

4.1. NGC 3227 

The original LAG estimate of the SMBH mass 
in NGC 3227 was ~2xl0 8 M Q (Salamanca et al. 
1994). Our value of (3.6±1.4) x 10 7 M Q is consis- 
tent with the mass of 4.9+^ x 10 7 M Q tabulated by 
WPM, which was based on an independent mon- 
itoring campaign (Winge et al. 1995). The value 
derived by Ho (1999) from H/3 was 3.8xl0 7 M Q , 
also consistent with our new SMBH mass in this 
galaxy. Schinnerer, Eckart, & Tacconi (2000) used 
CO gas kinematics to calculate a lower mass limit 
of ~T.5xl0 7 M Q enclosed within the central 25 pc. 
Again, this value is fully consistent with our re- 
sults. 

4.2. NGC 3516 

Kazanas & Nayakshin (2001) have noted the 
difficulty in reconciling the SMBH mass estimates 
for NGC 3516 arising from a variety of optical, 
UV, and X-ray constraints within certain models 
of AGN structure. Although our final estimate 
of M rev = (1.68±O.33)xlO 7 M lies slightly above 
a recent mass upper limit of (1.12±0.05) x 10 7 M Q 
derived from models of the X-ray warm absorber 
in this galaxy (Morales & Fabian 2002), other 
models of X-ray properties point towards a mass 
closer to 1O 8 M (Kazanas & Nayakshin 2001). 
Our results for this AGN are consistent with the 
value of ^2xl0 7 M Q from the earlier analysis of 
the LAG data by Wanders & Home (1994), and 



also matches well the mass of 2.3xl0 7 M Q Ho 
(1999) calculated using H/3 alone. Recent at- 
tempts at accretion disk modeling from X-ray data 
have given values of 3.1xl0 7 M Q (Czerny et al. 
2001) and 2xlO 7 M (Chiang 2002), both in rough 
agreement with the present work. 

4.3. NGC 4593 

For NGC 4593, we calculate that M rev = 
(6.6±5.2) x 10 6 M Q . An earlier reverberation map- 
ping campaign derived a mass of 2.2±^xlO 6 M 
from Lyman a (Santos-Lleo et al. 1995), and 
Kollatschny & Dietrich (1997) calculated M rev 
= 7xl0 6 M Q from the LAG Ha data (when cor- 
rected for our conversion between AV and the ve- 
locity dispersion). Ho (1999) estimated a slightly 
larger Hf3 time delay and found a SMBH mass of 
8.1xl0 6 M Q . With such large uncertainties, these 
differences cannot be considered statistically sig- 
nificant. 



5. THE M^-cr* RELATIONSHIP 

Since the discovery of the exceptionally tight 
correlation between a galaxy's SMBH mass and 
the central value of the velocity dispersion of 
the galaxy's spheroid (i.e. M^-cr*; Ferrarese & 
Merritt 2000; Gebhardt et al. 2000a), a number 
of attempts have been made to reconcile differ- 
ences in the derived slope of this relationship (e.g., 
Trcmaine ct al. 2002) as well as to develop a phys- 
ical framework that will naturally produce such a 
correlation. In addition, Gebhardt et al. (2000b) 
and Ferrarese et al. (2001) have examined whether 
SMBH masses derived by reverberation mapping 
are consistent with the M^-cr* relationship. Be- 
cause bulge velocity dispersion measurements ex- 
ist for all three of our targets, we are able to en- 
large the AGN sample used to investigate M^-cr* . 

Nelson & Whittle (1995) measured bulge ve- 
locity dispersions for 85 galaxies, including two of 
our targets. NGC 3227 was found to have er* = 
128±13 km s-\ and NGC 4593 has er* = 124±29 
km s _1 . Di Nella et al. (1995) published a value 
of cr* = 124±5 km s" 1 for NGC 3516. These val- 
ues for a-*, in conjunction with our derived M rev 
data, are consistent with the current fits to M^-cr* 
(Figure 4). This excellent agreement between the 
AGN and quiescent galaxy Mbh-<J* relationships is 
additional evidence that reverberation masses are 
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not subject to large systematic errors. 

In addition, we note the relative concentrations 
of the AGN and quiescent galaxy data in Figure 
4 toward the two ends of the relationship. This 
is a reflection of two facts: (1) in order to mea- 
sure cr* in AGNs, the continuum emission must 
be weak enough so as not to wash out the stellar 
absorption lines, implying the accessibility of only 
the dimmer AGNs, and therefore only the smaller 
SMBH masses; (2) the quiescent galaxies need a 
relatively large SMBH mass to allow the probing 
of the gravitational sphere of influence, whether 
by stellar dynamics or gas kinematics. The com- 
bination of these requirements gives rise to the ob- 
served distribution, in spite of the probable span of 
both galaxy types over the entire range of SMBH 
masses. 

6. SUMMARY 

As many of the forefront studies of AGNs are 
examining correlations between SMBH masses and 
properties of the AGN host galaxies, a uniform 
analysis methodology is important for minimizing 
systematic offsets between datasets. Time lags, 
calculated from the published light curves, were 
used in conjunction with rms velocity widths for 
optical hydrogen emission lines to estimate the re- 
verberation masses for NGC 3227, NGC 3516, and 
NGC 4593. These masses, listed in Table 4, can 
now be added to the compilation of Kaspi et al. 
(2000), as they were calculated in the same man- 
ner. In addition, these masses, when combined 
with published data on the host galaxy spheroid 
velocity dispersions, strengthen the claim that re- 
verberation masses are accurate measures of the 
true SMBH mass. 
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Jet Propulsion Laboratory, California Institute 
of Technology, under contract with the National 
Aeronautics and Space Administration. 
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Velocity (km s _1 ) 

Fig. 1. — Mean and rms Ha spectra for NGC 3227 
{(a) and (b)), NGC 3516 {(c) and (d)), and 
NGC 4593 {(e) and (f)). The y-axis scales of 
Fx arc arbitrary and the continua have been sub- 
tracted. 




Velocity (km s _1 ) 

Fig. 2.— Mean and rms H/3 spectra for NGC 3227 
((a) and (b)), NGC 3516 ((c) and (d)), and 
NGC 4593 ((e) and (f)). The y-axis scales of 
Fx are arbitrary and the continua have been sub- 
tracted. 
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Fig. 3. — Cross-correlation results for Ha (left 
panels) and H/3 (right panels) with continuum 
auto-correlation functions (dotted lines; cross- 
correlation of the continuum with itself) for 
NGC 3227 ((a) and (b)), NGC 3516 ((c) and (d)), 
and NGC 4593 ((e) and (f)). 
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10" 




Fig. 4. — The relationship between black-hole 
mass and the central value of host-galaxy bulge ve- 
locity dispersion for quiescent and active galaxies. 
The filled circles represent quiescent galaxies (data 
kindly provided by L. Ferrarese) and the filled 
triangles are AGNs from Ferrarese et al. (2001). 
The open triangles are the three AGNs discussed 
in this paper. The dashed line is the best fit to 
the quiescent-galaxy data, Mbh oc er 4,58 (Ferrarese 
2002). 
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Table 1 
Sampling Statistics 



Sampling Interval (days) 



Galaxy 


Datasct 


Number 


Average 


Median 






NGC 3227 


Ha 


23 


6.9 


5.1 


0.017 


1.367±0.155 




H/3 


14 


10.6 


10.9 


0.133 


1.597±0.181 


NGC 3516 


Ha 


18 


8.9 


8.8 


0.129 


1.541±0.065 




H/3 


18 


9.0 


8.0 


0.112 


1.474±0.063 




cont. 


35 


4.5 


3.0 


0.280 


4.027±0.364 


NGC 4593 


Ha 


22 


7.4 


5.0 


0.136 


1.653±0.069 




H/3 


11 


15.8 


11.7 


0.185 


1.996±0.102 
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Table 2 
Cross-Correlation Results 



Galaxy 


Emission Line 


Continuum 


Tcent 


rest 
'cent 


Tpeak 


—rest 
peak 






(A) 


(days) 


(days) 


(days) 


(days) 



NGC 3227 


Ha 


(24% 


stellar) 


5020 


25 q+ 14 - 9 


25 8+ 14 - 8 


25.0™ 


24.9±i 8 6 7 




H/3 


(24% 


stellar) 




12 -°±9 6 1 8 


12.0+fi 7 


13.0™ 


13.0™ 


NGC 3516 


Ha 






composite 


13.2; 


f 5.9 

-2.6 


i3-i±|;l 




13.7+H 




H/3 








7.4^ 


-5.4 
-2.6 


7.3±I;| 


10.2+°? 


10 l+° :9 


NGC 4593 


Ha 






6310 


4 6+ 2 - 5 




i.7il;| 


1.711;? 




H/3 








o 1+7.6 




q fi+i-2 

y -°-12.6 
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Table 3 
Velocity Width Results 



Ha H/3 



Galaxy 


Dataset 


AV 


AV rest 


Ay 


AV rest 






(10 3 km s- 1 ) 


(10 3 km s- 1 ) 


(10 3 km s- 1 ) 


(10 3 km s- 1 ) 


NGC 3227 


rms 


3.10±0.15 


3.09±0.15 


4.36±1.32 


4.34±1.31 


NGC 3516 




3.11±0.04 


3.08±0.04 


3.14±0.15 


3.11±0.15 


NGC 4593 




3.10±1.09 


3.07±1.08 


4.42±0.95 


4.38±0.94 


NGC 3227 


mean 


2.62±0.22 


2.61±0.22 


1.90±0.57 


1.89±0.57 


NGC 3516 




3.51±0.10 


3.48±0.10 


4.09±0.77 


4.05±0.76 


NGC 4593 




3.27±0.66 


3.24±0.65 


3.31±0.93 


3.28±0.92 
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Table 4 
Reverberation Masses 



Galaxy 


Ha-based 
(1O 6 M ) 


H/3-based 
(1O 6 M ) 


Combined 
(1O 6 M ) 


NGC 3227 
NGC 3516 
NGC 4593 


36.1™ 

i8.2tl:5 

6.3±|;« 


qo i +76.3 
°°- 1 -32.1 

10.3±^ 

o 7 +21.3 
S- 1-14.8 


36±14 

16.8±3.3 

6.6±5.2 
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